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Abstract 

That static electric and magnetic fields can store momentum may be per- 
Q \ plexing, but is necessary to ensure total conservation of momentum. Simple 



situations in which such field momentum is transferred to nearby bodies and 
point charges have often been considered for pedagogical purposes, normally 
lyj . assuming vacuum surroundings. If dielectric media are involved, however, 

,^ \ the analysis becomes more delicate, not least since one encounters the elec- 

^1 tromagnetic energy-momentum problem in matter, the 'Abraham- Minkowski 

r-| ' enigma', of what the momentum is of a photon in matter. We analyze the 

, ^ ' momentum balance in three nontrivial examples obeying azimuthal symme- 

try, showing how the momentum conservation is satisfied as the magnetic 
sj^ , field decays and momentum is transferred to bodies present. In the last of 

(^ \ the examples, that of point charge outside a dielectric sphere in an infinite 

2 ' magnetic field, we find that not all of the field momentum is transferred to 

■>!::;j- ■ the nearby bodies; a part of the momentum appears to vanish as momentum 

t^ . flux towards infinity. We discuss this and other surprising observations which 

O \ can be attributed to the assumption of magnetic fields of infinite extent. We 

emphasize how formal arguments of conserved quantities cannot determine 
which energy-momentum tensor is more "correct" , and each of our conser- 
Lj vation checks may be performed equally well in the Minkowski or Abraham 

'Vl ■ framework. 

. .=?, 

1. Introduction 

It may appear surprising at first sight that there is an electromagnetic 
field momentum with density 

g = D X B (1) 
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for crossed electric and magnetic fields in a vacuum field even under static 
conditions. This somewhat striking property of classical electromagnetic the- 
ory has actually served as a popular demonstration example to show how the 
law of momentum conservation works for a closed system. 

For illustration, let consider the setup of Fig. [H where there are two 
infinitely long metallic concentric cylinders carrying opposite charge with 
nearly the same radius a, one cylinder fitting inside the other. The following 
argument is adopted from the exposition in Aharonov and Rohrlich's book 
[l|. Initially the cylinders are made rotate in opposite senses, with constant 
and opposite angular velocities. At a distance x = xo from the common 
axis of the cylinders there is placed a heavy particle with charge Q. Assume 
that the cylinders gradually stop rotating, because of friction between them. 
What happens to the momentum balance? 




Figure 1: Two metallic cylinders rotating in opposite senses. 



Initially, if xq >■ a the electric field component E^ within the cylinder 
volume is practically the same as on the z axis (SI units assumed) 

F = ^ ^° (9) 

47r£o(a;2 + z2)3/2- ^^) 

(We assume here that the metallic walls have no influence on the electric 



field on the inside. This is clearly an unphysical assumption, but it has no 
importance for the argument. As an alternative, one might assume that the 
magnetic field were produced by external sources not influencing the electric 
field at all.) While the cylinders rotate at constant velocity the magnetic 
field inside is constant, Bz = Bq. The y component of momentum is thus 
Qy = —EqE^Bq, which implies that that the initial total field momentum Gy 
is 

Assume now that at t = the cylinder movement begins to slow down, for 
example due to friction. The magnetic field decays and eventually vanishes 
when the cylinders come to a halt. During this decay period an azimuthal 
induced electric field component £'™duced _ —Bo?/{2r) is induced in the 
exterior region. The total impulse imparted to Q by the force Fy = QE"^'^^'^'^'^ 
at r = xo is thus 

The expressions ([3]) and (j4]) are the same, thus confirming conservation of 
the y component of total momentum. 

Variations of the above argument, as already mentioned, are often pre- 
sented in connection with the foundations of electrodynamics - cf., for in- 
stance, Ref. [2|. These considerations actually have a close bearing also on 
the Aharonov-Bohm effect and the concept of "hidden momentum" in elec- 
trodynamics (cf., for instance, Refs. |3[; a recent feature article is Ref. [4]). 

The question now emerges - and this is the topic to be considered in 
the present paper - in which ways does the argument above change if the 
test region is filled with a dielectric medium? As we will see, analysis of 
such situations can be quite subtle. That some complications necessarily 
have to occur, is evident already from our initial expression ([T]) for the mo- 
mentum density: in the presence of a medium that expression is no longer 
trivial but becomes the same as the Minkowski expression |5j], just one of 
the alternatives that is frequently discussed in the current literature. Other 
alternatives, such as the Abraham expression [6| which is also well known, 
would correspond to different expressions for g and thus modify the momen- 
tum balance calculation. The modification that we are considering, is thus 
closely connected with the problem what has turned out to be known as the 
'Abraham- Minkowski problem'. 



A nice introduction to the Abraham-Minkowski problem is found in M0ller's 
book 7|. A review article was given by one of the present authors some years 
ago i&\. Another review, up to 2007, is given by Pfeifer et al. [9!]. Some more 
recent papers are Refs. [lO|, lUl, Ijji ll3|, [ij, ll5|, ll6|, [22 . 



We will consider three different examples of momentum conservation 
checks in the presence of dielectric media, in order of increasing complex- 
ity. First, we consider a natural extension of the introductory example given 
above, that of a point charge in the vicinity of a magnetic field which is 
restricted to a cylindrical volume, all embedded in a dielectric medium of 
spatially constant permittivity. 

In the light of this first example we discuss an important point, namely 
that, contrary to the claims by some authors in the past, such formal consider- 
ations of conserved quantities cannot be used to determine the "correctness" 
of one formulation of electromagnetic momentum over another, for example 
the famous formulations due to Abraham and Minkowski mentioned above. 
Next, in Sect. 3 we consider the case of an infinite dielectric cylinder with 
a parallel line charge placed on the outside. The final, and most involved, 
example, is the geometry of point charge in the vicinity of a dielectric sphere 
placed in a magnetic field of infinite extent. In this final example we ob- 
serve that not all of the field momentum is in fact transferred to the sphere 
and point charge as the magnetic field decays: a part of it is transmitted, 
at least formally, as momentum flux towards infinity. The point is a subtle 
one, however, and paradoxes arise due to the unphysical assumption that the 
magnetic field fills all of space, which we discuss before concluding. 

2. Example I: Point charge near confined magnetic field in a medium 

Consider the situation sketched in Fig. |2) a point particle of charge Q is 
placed in the exterior vicinity of a long cylindrical volume of radius h wherein 
a uniform and axially directed magnetic field B = BqTj is imposed. At a point 
in time the magnetic field starts decaying, and we study how the momentum 
of the electric and magnetic fields in the cylindrical volume is transferred 
to the particle. This is essentially the same system as in the introductory 
example, but in the presence of an infinite medium. We assume there is 
zero magnetic field at radii greater than h. For generality let the system be 
embedded in an infinite medium of permittivity e (we define the permittivity 
e as nondimensional, so that the static constitutive relation reads D = esoE). 



To facilitate calculations, let us place the system in a large perfectly 
conducting cylindrical "box" of length L and radius R. We will assume 
eventually that L and R are much greater than other length scales in the 
geometry. The electric potential at the conducting surfaces aX r = R and 
z = ±L/2 is $ = 0. The charge Q sits on the x axis at position r' = vq = 
Xq, 9' = 0, z' = 0. We initially assume the particle to lie outside the area of 
the magnetic field, vq > b. The introduction of the large external cylindrical 
"box" enables us to deal with a completely closed physical system. 




Figure 2: Charge near magnetic field confined to cylindrical volume, all embedded in an 
infinite medium of permittivity s. 



One may start the calculation by use of the Green function G for the 
Poisson equation, which satisfies 



V'G'(r, r') = -—6{r - r')6{e - e')6{z - z'). 
r 



If G is known, the potential follows from 



G. 



(5) 



(6) 
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The solution of the potential problem in the case of a vacuum filled closed 
cylinder of length L and radius R can be found, for instance, in Ref. 17 



When generalized to a medium, the solution can be written 
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Here /„ and K^ are modified Bessel functions, r> (r<) is the greater (smaller) 
of r and tq, and we have defined q as 
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We require the r and 6' components of the electric field E = — V$ for radii 
< b (i.e. r < ro): 
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A prime means differentiation with respect to the whole argument. 

We shall calculate the y component of the total Minkowski momentum, 
similarly as above (from the symmetry it is obvious that the electromagnetic 
field can have no net momentum in the x or z directions). This momentum 
is 

G^ = Bo rdr dz d9i-DrCos9 + DesinO). (11) 

Jo Jo Jo 

The integrals may be processed by standard relationships including 



/ d9cos9cosm9= / d9msm9 sin 1719 = iridmi + 5m,-i), 
Jo Jo 
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and 



/ drii (pqr) = — [Iq (pqb) - 1] ; 
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noting that Ii{x) = /_i(x) and Ki{x) = K_i{x). After some algebra we 
obtain the expression, inserting rg = Xq, 



G = > - sm'^ — 1 1 [pqb) Ki [pqxo 

p=l ^ 
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This expression gets contributions from p = 1,3,5.... It is therefore conve- 
nient to introduce a new integer I = {p+l)/2 such that all values / 
contribute. In terms of / and q = tt/L we then get 
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This expression holds for all values of 6 up to 6 = xq. 

We now take the limits R,L ^ XQ,b. Whichever way this combined 
limit is taken, one finds that the second term of the square brackets in (1T6|) 
vanishes. Moreover, by using the expressions Ii{x) ~ x/2, Ki{x) ~ 1/x 
which hold for small arguments x, we get Ii[{2l — l)qb]Ki[{2l — l)qxo] ~ 
b/{2xo) for moderate values of I, which is in order since the resulting / series 



is convergent. Observing that 1 — | + | 
limit the very simple expression 
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+ ... = 7r/4 we obtain in this 



(17) 



The result is independent of the large external cylindrical box as it should. 
Moreover it is identical to Eq. (jl]), although we have not made any assumption 
here that Xq be much greater than b. 

Consider now the time interval t > in which the magnetic field decays. 
Since the outer cylinder is far away we can assume that the induced electric 



pinduced __ ^ rj 



field component E^^'^^'^°'^ at radius tq is not influenced signiflcantly by the 
condition that it must vanish on the boundary at r = i? ^ rg. The induced 
azimuthal electric fleld ii;™<i^ced g^^. ^ radius h < r <^ Ris then simply found 
from Faraday's law as 

The y component of the impulse imparted to Q because of £;^°duced jg^ ^^ -^^ 
Eq. (jl]), the time integral of Fy = QEg evaluated at r = xq, 9 = 0, 

Equations (TT7|) and (TT9|) are in agreement, thus demonstrating the momentum 
balance. 

We have restricted the width of the uniform initial magnetic field so that 
b < Xo- We will consider the case b > xq below in section 12.21 but let us 
first turn to an important point regarding different proposed forms of the 
electromagnetic momentum in media, which we are now able to illustrate in 
the light of the above example. 

2.1. Abraham vs. Minkowski momentum 

We have so far considered the Minkowski momentum only. As mentioned 
in the Introduction there exist other proposals for the momentum also. We 
will argue here how, in contradistinction to some claims in the literature, one 
cannot use momentum balance considerations such as those we consider in 
this paper to determine which energy-momentum tensor is "correct" . 

The most well-known among the alternatives to Minkowski's tensor is 
that introduced by Abraham [6|. In the Abraham tensor, the momentum 
density g^ is given as 



1 



r^ = ^(E X H). (20) 



Equations flT6|) or fITTj) would in this case be replaced by 

1 



G^ = -G^\ (21) 
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The simple relationship between expressions flT7|) and (1191) is thus violated, 
Gy ^ Py, since Py = G^ . The question naturally emerges: does this imply 
that the Abraham form runs into conflict with the conservation of total 



momentum? The answer is no. The reason is that the Abraham energy- 
momentum tensor describes a non-closed system and is not divergence- free. 
There exists an Abraham force density f"^, acting on the medium even in 
spatially uniform regions. Its magnitude is 

f^ = ^|(E X H) (22) 

It is usually called the Abraham term. Under static conditions corresponding 
to t < the term plays no role at all, but for t > when B 7^ the force 
(!22|) has to impart a mechanical impulse to the medium. The impulse thus 
imparted is (we assume non-magnetic medium throughout for simplicity) 

AG^ = f dt f rfVf^ = -i^ f rfV(Do X Bo) 
Jo Jv ^ Jv 

=G^ - G^ (23) 

where V is the volume in which both the electric and magnetic fields are 
nonzero, and subscript means the value at t = 0. 

As Eq. ( I23l) shows, the difference between the Abraham and Minkowski 
values of the initial field momentum is imparted to the integration volume 
during the decay period, and thus the momentum balance is ensured also in 
this formalism: the total momentum transferred from the initial fields to all 
bodies surrounding bodies is Py + AG^ which equals the initial momentum 
in both formalisms. Thus the momentum conservation is demonstrated also 
in the Abraham framework. 

There exist several papers in the literature making unjustifiably strong 
conclusions based upon analyses of conserved quantities, such as total mo- 
mentum. One example is that of Skobel'tsyn pLSj] (that paper was discussed 
also in the review [Sj, page 192). There is generally no definite answer about 
"correctness" to be obtained from this kind of analysis; one ends up perform- 
ing consistency checks of the formalism rather than giving real derivations. 
In our opinion it is instead from experimental information that one can get 
reliable information about appropriateness and usefulness (rather than "cor- 
rectness"), of the various alternatives for electromagnetic momentum. 

2.2. Second case: b > xq 

Let us briefly return to Example I, but now allow the charge to lie 
within the cylindrical region in which the ^-directed magnetic field is initially 

9 



nonzero. Since the calculations are closely analogous to those in section |21 
we do not here give the full details of calculations. The contribution from 
the region r < xq is 

Gf - = ^QBoXo (24) 

as seen immediately from Eq. ( TTTl) . To calculate the momentum contained 
within the annular cylindrical region xq < t < 6, we require the electric 
fields in this region, easily obtained by differentiation of ([7]). Performing the 
integrals over z and 6, and summing over m (only m = ±1 contribute as 
before) we obtain the annulus momentum 



or =^ E-^' (?) ^i(P^-o) / rdr L 

p—l J Xq K. 



1 
+ — 

qr 



.^^^^^^^^-f^'^^^^^). 



Kiipqr) - -— — —I^ipqr) 
h{pqR) 

(25) 



Solving the integral over r as above and taking the limit L,R ^ h, Xq we 
obtain 

Gf+= hm ^^5M^^^v[;(i?/L) ^^°° 0. (26) 



[The function ^(x) = E/=i(-l)'(2^ - ^fKi[{2l - l)x]//i[(2/ - l)x] exists for 
X >0]. 

We thus see that there is no contribution of momentum from the annulus 
of radii > xq and the total momentum is given by flM|) . Note that this also 
holds in the limit h — )■ oo, a point we will discuss later. One sees almost 
immediately that the total momentum is again transferred to the line charge 
and is conserved. 

3. Example II: Infinite dielectric cylinder, with external parallel 
line charge 

Consider as next example a dielectric nonmagnetic cylinder of radius a, of 
infinite length and permittivity e, as shown in Fig. 121 We employ cylindrical 
coordinates r, 9, z, where the z axis coincides with the cylinder's symmetry 
axis. In the outside region r > a we assume that there is an infinite parallel 
line, of charge A per unit length, parallel to the cylinder. The line passes 
through the xy plane at position x = Xq, 9 = 0. The electric field lies 
entirely in the xy plane, E = E^r + EqO. For times t < there is a uniform 

10 



magnetic field Bq z imposed in the z direction. We will assume as in Example 
I that the magnetic field is restricted to radii smaller than a radius b. The 
three cases b < a and a < b < xq and b > xq will be treated separately. As 
before, we assume that the magnetic field decays for t > 0. 

We shall examine the momentum balance in the transverse y direction, 
per unit length of the cylinder. As the length is infinite, this means that we 
are effectively dealing with a closed physical system. 




Figure 3: A line charge A per unit length outside an infinitely long dielectric cylinder 
of radius a, positioned coaxially with a magnetic field confined at radii < b. The case 
a < b < Xq (second case) is shown here. 



3.1. First case: b < a 

Consider the internal electric field E^. The potential $~ is the same as 
that produced in a homogeneous medium (with permittivity e) by a fictitious 
charge 2Xe/{e + 1) placed at Xq |19| : 



$ = 



In Tr^ — 2rxr) cos 9 + x^] . 

2-Keo e + l ^ " °J 



(27) 
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It leads to 



A 1 r — xq cos 9 

tieq £ + 1 r^ — 2rxo cos 9 + x^' 

A 1 Xq sin 9 

neo e + 1 r^ — 2rxo cos 9 + Xq 



E; -r^-^ „, Z7:Z:. . . (28a) 



Ei =^ T^ ,., :r:..^ - (^ 



The 2/ component of the interior momentum density is 



from which it follows that 



g^- = -Bo{D; cos 9- D, sin^), (29) 



Gf'= r d9 frdrgf-. (30) 

Jo Jo 



We get 



_ AB,^ r- t l-(fo);°^f ,31) 



TTXo £ + 1 7o 7o 1 - 2(r/a;o) cos6' + (r/xo)^ 

Recognising the Fourier series 



1 — s cos 9 
1 -2scos^ + s2 



^s"cosn^, s < 1, (32) 



n=0 



we see that the integral over 9 in Eq. fj3T|) gets contribution only from n = 0. 
Thus 

One may note that in the limit e ^ 1 this expression becomes identical to 
dlj) and (fTTll . with A replacing the point charge Q. 

Consider now the magnetic field decay period t > 0, in which we have 
to include the induced electric field ii;™duccd _ —i)'^j^^2r for r > b. There is 
no electromagnetic force acting in the homogeneous interior; all forces are 
surface forces. Such forces can be calculated by taking the difference between 
the radial components of the Maxwell stress tensor on the outside and the 
inside. However, it is in our case more convenient to start from the fact that 
the force density in the surface layer is —{eo/2)E^'Ve. By integrating this 
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across the surface layer and using electromagnetic boundary conditions we 
can express the surface force density cr as 

<"r, d) = \eo{e - 1) [E'^ir, 9) + eE',{r, 9)] r, (34) 

Z a— 

thus showing that we need the interior fields only. The transverse y compo- 
nent is 

^2/W = ^eo{e-l)[E^g{a',9)+eE^^{a-,9)]sm9, (35) 

where the internal fields are given by (128|) with the addition of _E'™duced ^^ 
the 9 component. When inserting these fields into the expression for the y 
component of force on the cylinder, 

^cylinder _^/ d9ay{a,9), (36) 

Jo 

we do not need to evaluate all the individual terms since only those containing 
B can contribute. Indeed the only contribution comes from the cross term 
in {EqY, and we obtain 

^cylinder \BP S - 1 f'^ d9 sin' 9 



^ 2ttxo e + I Jq 1 — 2{a/xo) cos9 + (a/xo)" 

A5fe2 e - 1 



2x0 £ + 1 
where we use the Fourier series 



(37) 



oo 



!!!i4 o = > ^ s"-' sin n9 (38) 

l-2scos^ + s2 ^ ^ ' 

n=l 

for s < 1, which upon insertion into our integral contributes only for n = 1. 
This leads to the momentum transfer 

pcylinder / ^j. ^cylinder AIJQi)t ± (QQ\ 

In addition to the force on the cylinder surface there is also an azimuthal 

y 

ippH fiplrl i« ncrnm Rinduced ^ _ 

2r 



13 



force Fy"^^ = A-E™'^"'^'"^ acting on the line charge A during the decay period. 
The induced field is again i^;™duced _ — |-ij go, with r = xq, 

Fr=-^, (40) 



leading to the impulse transfer 



iline 



dtF 



line 



2xn 



(41) 



The total impulse transfer to the system during the decay period is accord- 
ingly 

p _ pcylinder , pline _ ^-°0" ^ (An\ 

^'~^y ^^y ~ xo e + i ^^^' 

This expression is seen to be the same as Eq. f l33p . The momentum balance is 
thus verified, within the framework of the Minkowski formulation. To check 
the momentum balance in the Abraham case, one would have to take into 
account the Abraham force acting on the interior of the medium during the 
decay period, similarly as we did above. 

3.2. Second case: a < b < Xq 

For the case 6 > a we must also consider the electric field from the line 
source in the region r > a. The potential $^ is that produced in a vacuum by 
the actual charge A, plus that of two fictitious line charges, one of magnitude 
A' = — A(e — l)/{e + 1) located at the position xi = o^/xq, the other of 
magnitude —A' located at the origin 19|, wherewith: 



$" 



A 



iner 



ln(r^ — 2rxo cos 9 + Xq) ln(r^ — 2rxi cos 9 + xf 



e + 1 



+ 



Inr 



e + 1 
This gives the fields 



(43) 



e: 



Ei 



A 



27reo 
A 



r — Xq cos 9 



— 2rxo cos 9 + Xq 



Xq sin 9 



e-1 
e+1 
€-1 



r — xi cos 9 



- 2rxi cos 9 + xf 
Xi sm9 



r^ — 2rxQ cos 9 + x^ e + Ir'^ — 2rxi cos 9 + x\ 



The resulting momentum density is found from (1291) to be 



9^^ 



r,9) 



ASo 

2tt 



rcos9 



Xq 



e-1 



r cos 9 — xi 



2rxo cos 9 + Xq 



£ + Ir"^ — 2rxi cos 9 + x\ 
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e — 1 cos 9 

e + 1 r 



(44) 
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where /3i = h/xi and fi = r/h. The contribution to G^^^ is now 

nb r27r 

0^'+= / rdr / d9g^'+{r,9). (45) 

The first term of (jUj) may be integrated over 6 as before, using f p2|) . Since 
Xi < r in the annular region between radii a and 6, however, we must use a 
different series expansion for the second term. 



COS y - s _ ^ ^2n-i ^^g 2n^ + COS ^ 5^(-l)"s'" 



2^cos2A;e 



fc=i 



l-2scose + s2 

n=l n=0 

= COS 6* + s COS 26^ - s^cos6'(l - 2cos26') + s^cos46' + ... 

(46) 

for s < 1, from which we see that the second term of (jUj) gives no contribu- 
tion when integrated from ^ = to 27r. The same is clearly the case for the 
last term in Eq. (jS]), so that only the first term of (jH]) contributes; 

Gf + = /' rdr r dOg^^r, 6) = ^^°^f ~ ^'^ (47) 

Ja Jo ^^0 

After inserting the momentum from the internal region from ( l33l) with b ^ a, 
we find the total field momentum to be in the Minkowski formalism 



G,^Gr^Gr-j^[^'^j^a'). (48) 

We again turn to the period in which the magnetic field decays. The 
momentum transfer to the line charge and cylinder, respectively, are simply 
evaluated just like before, to yield 



2a;o 



pline ^0_. (49) 



pcylinder _ A-Dpa^ £ 1 , , 

Adding Py^'^ and p^yimder ^g retrieve the expression f HHj) . and the momentum 
conservation is once again ascertained. 
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3.3. Third case: b > Xq 

Consider finally the case where the line source lies within the volume 
in which the magnetic field is present. The momentum density inside the 
sphere is again given by fl55]) with 6 — > a, and the momentum density for 
r > a again given by (jS]). As in the previous case, only the first term of 
(jH]) contributes. Similar to the calculation in the previous section, different 
Taylor expansions of the fraction (r cos 6 — xo){r'^ — 2rxo cos 6 + Xq) must be 
employed depending on whether r < Xq or the other way around: the series' 
( 15^ and (H^ . respectively. The latter of these, as previously argued, gives 
no contribution after integration with respect to 6, so we conclude that the 
annular region xq < r < b, just as in the second case of Example I, carries 
no momentum, and the total field momentum for t < is found as 

One easily verifies that the to terms in the paranthesis exactly equal the mo- 
mentum transferred to the line source and dielectric cylinder, respectively. 
Note once again that this result holds also in the limit where 6 — )■ oo. Sub- 
tleties connected to this limit are discussed in section 

4. Example III: Dielectric sphere in an infinite magnetic field 

In the above two examples, we had the opportunity to exploit the cylinder 
symmetry property. As a final example we shall consider briefiy also the 
analogous problem with a dielectric sphere of radius a, acted upon by an 
external charge Q, and also exposed to a uniform magnetic field which for 
reasons of tractability we shall assume to be of infinite extent. 

Employing usual spherical coordinates r, 6, cj) it is mathematically conve- 
nient to place the external charge Q on the z axis, at position zq> a and let 
the transverse uniform magnetic field be directed normal to the sphere-charge 
line, along the x axis: Bq x for times t < 0. As mentioned, the magnetic field 
is present everywhere, on the inside as well as on the outside. 

From Ref. [20[ we have for the potential on the inside 
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Figure 4: A point charge outside a dielectric sphere in an infinite x-directed magnetic 
field. 



giving rise to the field components 



e: 



Q ^ 21 + 1 Ir^-' ^, 



'0 ~ '-t- + ' + 1 Zq 



E- =-^ y ^i±l^ 4tt p/(cos ^) sin ^, 

E-^ =0. 
The y component of total momentum in the sphere is 



(53) 

(54) 
(55) 



G'f- = 27r / r^dr / sinOdOg^ 



y ' 



where the density for t < is 



Qy^ = Bq{D^. cos 9 — Dq sin I 



Inserting the expressions we get 



G 



M- 



>"^^t {i) 



l+l 



2/ + 1 



(/ + 2)(fe + / + l) 



X / [/cos0Pz(cos^)-sin^P/(cose)]sine(ie. 
Jo 



(56) 



(57) 



(58) 
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, 2 

a \ e 



Some algebra leads to the expression 

Gf- ^ -QB,a ^^j ^, (69) 

We have to consider the y component of momentum from the fields on 



the outside as well. The potential $"•" in that region is 20 1 



which gives for the fields 

Q r-zocos9 Q >^ /(/ + !) a'^+^ PijcosO) 



47r£:o ^^ 47reo T^f ^^ + ^ + ^ 4"^^ 



EJ =0, (63) 

where ^ = (r^ — 2rzocos^ + 2:0)^^^ i^ ^^^ distance between Q and the field 
point. These expressions are inserted into the analog of Eq. ( 1571) in the outer 
space (superscript minus replaced by plus). An explicit calculation shows 
that the contribution from the medium-generated fields, i.e. the last terms 
in Eqs. ( 16T|) and ( l62l) . is actually zero. The contribution from the first terms, 
which stem from the point charge alone, are calculated by means of the 

relation 

sin6'(r cos6' — 2:0) 2^, , , . 

d9 . , , 2 , -^ = — 26(^0 - r) (64) 



+ Zq — 2rzo cos 6*)^/^ z. 



where Q{x) is the unit step function. We obtain the result 



or = -^{z-^ - a-^) (65) 



so the total momentum of the electromagnetic fields before the decay period 
is 

A note is warranted at this point. The second term of fl66|) is in a sense 
arbitrary: it refers to the absolute position of the point charge only, and not 



the sphere. The first term, in comparison, depends only on parameters which 
are independent of the coordinate system: the sphere radius and the ratio of 
the radius to the distance of the point charge from the sphere center. Would 
not then the second term change its value upon a linear transformation of 
coordinates, but not the first? And secondly: how can a term not referring to 
the sphere possibly be non-zero, when the remaining system of point charge 
and infinite field has perfect spherical symmetry? 

In fact we have seen this paradox also in the previous chapters, although 
it was not discussed at the time. Both the field momentum (124|) for the 
particle in a cylindrical field and (I5T]) are valid in the limit where the magnetic 
field fills all of space, although were one to take that limit before doing the 
calculations, one would get the result zero because of symmetry. Also here 
we must interpret the infinite magnetic field to be the limit where some 
confinement volume is taken to infinity. Depending on the original symmetry, 
the charge alone can constitute a symmetry break (provided the limit of 
infinite field is taken at the end) and give rise to a field momentum different 
from zero, although the physical meaning of this momentum is not so obvious. 
We return to this point in Section O 

Consider now the magnetic decay period, t > 0. The impulse trans- 
ferred to the sphere, calculable from the fields at r = a—, necessarily has to 
correspond to the surface force density component 



(^y = 7^eo{e~l)[El + eE', 



sin 6 sin i 



(67) 



As discussed we assume azimuthal symmetry for the induced electric field 
■gjmducod g^j.Q^]2d the direction of B. One finds 



E 



induced 



-rB'k X f, 



(68) 



where f is the unit vector directed normal to, and away from, the x axis 
[here f = z/ sin a, and f = cos ay + sin az wherein cot a = y/z = tan 9 sin 0] , 
implying the components 



pinaucea ri 



linduced 
T-iinduced 



E. 



induced 



2 
1 



-rB sin i 



-rB cos 9 cos ( 



(69) 
(70) 

(71) 
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Letting r = a and adding these expressions to Eqs. (135]) - fl55p for the interior 
fields at r = a—, we obtain after some algebra the following impulse transfer 
to the sphere 

Pf"' = -\QB,a{£J'-ll. (72) 

The momentum transferred to the point charge is also easily calculated, 
assuming the charge stays close to 6* = throughout the decay period: 

p^charge ^ J_QB,z,. (73) 

This is the same result obtained in Example I, case 2 above, Eq. f ll7p . and 
for obvious reasons: in both cases the induced electric field is calculated from 
the single property of cylindrical symmetry of B around the x axis. 

There appears to be a discrepancy present: unlike in the previous ex- 
amples, the momentum transferred to charge and body does not equal that 
initially present in the electromagnetic fields. There is a "missing" momen- 
tum 

^Gy = --QB^a{-] i^—+QB^z,. (74) 

D \Zq J e + Z D 

Where has the momentum gone? The only possibility is that it has escaped 
the system towards infinity. We will now verify that this is indeed so. 

This we do by calculating the net fiux of y-directed momentum through 
a control surface which we place far from the sphere and point charge. We 
let the control surface be spherical for simplicity, with radius i? 3> a, zq. The 
flux of y-momentum through the surface is given as 

Gy = -i y-T-dA = -R^ desinO d(j){y -T ■ r){e,(j)) (75) 
Jcs Jo Jo 

where the stress tensor in vacuum is 

T = eoE®E+— B(8)B-^(E2 + c^B2)I (76) 

/io 2 

where I is the unit matrix (Kronecker's delta). In vacuum the stress tensor 
is related to the (any) energy momentum 4— tensor S^u by Tik = —Sik, latin 
letters denoting spatial dimensions. We obtain 

y • T • f = — [sin 6 sin (f)E'^ + 2 cos 9 sin (pEgEj. — sin sin (jyE"^ 

+ cos (pE^Ej. — sin 6 sin cpE"^ — sin 6 sin (j)(?Bl] . (77) 
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Only a few of these terms give a contribution after integration with respect 
to 0. The field expressions are simplified after expanding in powers of a/ R 
and zq/R 



r-2-K 

/ t/0(y • T • f ) 



QBR 



2cos^ Zq , . 

+ ?r^(3 + 5cos20) 



£ + 2 zlR^ 



i?2 2i?3' 

3 

(4 cos^ Q - sin^ 



'^ '■ .2/1 „• 2, 



/ I 



-n^4i- (^«' 



The higher order terms will vanish when taking the radius R to infinity and 
we need not consider them. Finally, performing the Q integral we obtain the 
momentum which has disappeared as flux during the decay period 



oo 



ar^l .*G. ^ f^ (.„ - i^l) (79) 

(recall: Jq dtB = —Bq) which implies that exactly all of the "missing mo- 
mentum" ^Gy in fl71|) has disappeared to infinity. In other words, a momen- 
tum has been imparted on the infinite background medium during the decay 
period, not dissimilar to the effect of the Abraham force in section 12. 1[ This 
striking result is discussed in the following. 

5. Discussion and conclusions 

We have verified the momentum balance of electromagnetic momentum 
in decaying magnetic fields in the presence of dielectric media in three dif- 
ferent configurations: first a point charge and a homogeneous magnetic field, 
the latter restricted to a cylindrical volume, embedded in an infinite dielec- 
tric medium, secondly a point charge in vacuum outside a dielectric cylinder 
aligned coaxially with the same magnetic field, and finally the case of a dielec- 
tric sphere and a point charge in the presence of an infinite and homogeneous 
magnetic field. While the momentum balance is simply and uniquely calcu- 
lated in the first two examples, the third example is more subtle, and in that 
case we find that some of the momentum disappears to infinity as momentum 
flux, calculated by use of the energy-momentum tensor. 

The fact that part of the initial momentum leaves the system is a striking 
find. The "physicalness" of the result is not so obvious, however, because the 
infinite extent of the magnetic field introduces arbitrariness, as we will now 
explain and discuss. First, however, note a second and similar paradox which 
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appears upon comparing Examples I and III: In Eq. fl66p we found one term 
which clearly refers to both sphere and charge, and another referring to the 
charge alone. It is natural, as discussed above, to interpret the second term 
as pertaining to the single charge and the (axisymmetric) magnetic field, and 
would if so be expected to have the same value even if the sphere were not 
present. Surprisingly, however, the latter term does not equal that of just a 
charge inside a cylindrical magnetic field, Eq. fl2^ . even when the cylindrical 
field volume radius is taken to infinity. 

Both of these surprising observations are connected to the unphysical 
assumptions of a magnetic field which fills all of space, which we made in 
order to be able to perform explicit calculations. In fact, this limit is not 
unique - it matters what the shape of the field volume is as the limit is taken. 
Note how, in this limit h — > oo, the results Eq. f l2^ and the first term of f lSTj) 
(just like the second term of Eq. (|66|) ) become arbitrary: once the magnetic 
field is infinite in every direction, new symmetries emerge. We are now free 
to move and rotate the coordinate system freely. Then the formal momentum 
of the electromagnetic field of a point charge in an infinite magnetic field, 
which depends on the absolute coordinates of the charge, can take any real 
value, depending on how the limit is taken! 




(b) 







Figure 5: Gedanken experiment: two different "field volumes" in an infinite medium. 



To illustrate even further, consider the following gedanken experiment, 
illustrated in Fig O The first situation. Fig. \^, is the same as in Example 
I, the second. Fig. [5]d, is the hypothetical situation of a magnetic field Bqz 
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confined to a spherical volume of radius b. The latter situation is clearly 
unphysical since magnetic field lines must be closed loops, yet let us ignore 
this fact for now. Using the potential of Eq. fl60|) we can calcinate (with a 
cyclic permutation of coordinates) the initial momentum contained in the 
spherical volume. We get 

G^ = X < ^0 ' 'for the cylinder, and 

2 \ zo, b> zo 



tM 



QBq J ^, 6 < zo, 



G,. = X < ^0 for the sphere. 

^ ^ ^ Zo, b>Zo 

The former is the result of Example I, the latter we recognize as the momen- 
tum term in the spherical example, Eq. fl66|) . which refers only to the charge. 
When 6 ^ oo in both cases one obtains the magnetic field filling all of space, 
but formally the momentum arising from the point charge is different in the 
two cases. 

The discrepancy between the initial momentum and that transferred to 
the sphere and point charge in example III then has the following explana- 
tion: the former is calculated taking the infinite magnetic field to be the limit 
of the spherical field in Fig. 5, the latter taking the cylindrical limit. The 
difference, being purely a mathematical artifact, cannot give rise to physi- 
cally measurable forces and is therefore absorbed by the infinite background 
medium. 

The subtlety relating to infinite magnetic field extent can be tackled by 
prescription: One may redefine the initial momentum of the initial fields by 
subtracting the "unobservable" amount which disappears as flux when the 
magnetic field decays (calculating both quantities in the same formalism). In 
Example III we then get the "renormalized" momentum ~ 22° fr^ ~ ^°^° , 
which has the "cylinder" prefactor 1/2, which is satisfactory since the cylin- 
drical limit in Fig. |5]is the more physical (such a field can be produced inside 
a large solenoid, the magnetic field lines closed at infinity). This "renormal- 
ization" procedure is closely analogous to the transfer between Minkowski 
versus Abraham momentum which we discussed in section 12. It there the 
Abraham momentum is obtained by subtracting from the Minkowski mo- 
mentum the impulse transferred by the Abraham force during the decay 
period. It is tempting to state that the "renormalized" momentum is the 
more physical, although as the century-old Abraham-Minkowski controversy 
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has shown that such statements are subtle. The momentum balance in the 
presence of an infinite magnetic field, and hence the formal consistency, is 
found to be in order with or without such "renormalization" as it should. 

We emphasize once again that, while we have used the Minkowski for- 
malism for electromagnetic momentum in matter in the examples calcu- 
ated in previous sections, the Abraham momentum could equally well have 
been made use of as discussed in section 12. 1[ Considerations of conserved 
quantities cannot be used to support a particular energy-momentum tensor 
over another: such calculations are merely internal consistency checks which 
both the Abraham, Minkowski and other proposed energy-momentum ten- 
sors should all pass if calculations are carried out correctly. It is our belief 
that the appropriate criterion for distinguishing between tensors is whether 
they can explain experimental observations in a straightforward manner. A 



classic example of such an experiment is that of Walker and Lahoz [21 



Another set-up was recently suggested by us [22 
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